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Abstract. Many clues about the galaxy assembly process lurk in the faint outer regions of 
galaxies. Although quantitative study of these parts has been severely limited in the past, 
breakthroughs are now being made thanks to the combination of wide-area star counts, deep HST 
imagery and 8-m class spectroscopy. I highlight here some recent progress made on deciphering 
the fossil record encoded in the outskirts of our nearest large neighbours, M31 and M33. 
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1. The Faint Outer Regions of Galaxies: Motivation and History 

The study of the faint outskirts of galaxies has become increasingly important in re- 
cent years. From a theoretical perspective, it has been realised that many important 
clues about the galaxy assembly process should lie buried in these parts. Cosmological 
simulations of disk galaxy formation incorporating baryons have been carried out by sev- 
eral groups and yield predictions for the large-scale structure and stellar content at large 
radii - for example, the abundance and nature of stellar substructure and the ubiquity, 
structure and content of stellar halos and thick disks. These models generally predict a 
\q ' wealth of (sub)structure at levels of jiy ~ 30 mag/D" and fainter; their verification thus 
requires imagery and spectroscopy of galaxies to ultra-faint surface brightness levels. 

Since Malin first applied his photographic stacking and amplification technique (e.g. 
Malin et al. 1983), it has been known that some galaxies possess unusual low surface 
brightness (LSB) structures - shells, loops, asymmetric envelopes - in their outer regions. 
Although limited to /is < 28 mag/D", these images were sufficient to demonstrate that 
even the most "normal" nearby galaxies could become very abnormal when viewed at 
faint light levels (e.g. Weil et al 1997). Follow-up study of Malin's LSB features, and 
several more recently-discovered examples (e.g. Shang et al. 1999), has been severely 
• <-h , limited due to the technical difficulties associated with detecting and quantifying diffuse 
light at surface brightnesses ~ 10 magnitudes below sky. The currently most viable 
technique to probe the very low surface brightness regions of galaxies is that of resolved 
star counts (e.g. Pritchet & van den Bergh 1994); I review here some recent results from 
studies of the faint outer regions of M31 and M33. 



2. Substructure Around Galaxies 

Figure 1 shows maps of the red giant branch (RGB) population (age > 1 — 2 Gyr) in 
M31 and M33, covering 40 and 7 square degrees respectively. Although the maps reach 
the same limiting depth (w 3 magnitudes below the tip of the RGB), it can be clearly 
seen that outer regions of M31 contain a wealth of substructure that is not present in 
M33 (Ferguson et al. 2002, Ferguson et al. 2007, in prep). The faintest features visible by 
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Figure 1. INT/WFC RGB star count maps of M31 (left) and M33 (right). These maps span 
125 x 95 kpc 2 and 56 x 56 kpc 2 on a side respectively. 

eye in the M31 map have effective V-band surface brightnesses of [iy ~ 28 — 31 mag/D"; 
if substructure is present in M33, it must be of lower surface brightness than this. 

Various features can be seen around M31, including a giant stream in the south- 
east, stellar overdensities at large radii along the major axis, a diffuse structure in the 
north-east and a loop of stars projected near NGC 205. It is of obvious importance to 
establish the origin of this substructure. Bullock & Johnston (2005) predict substructure 
in the outer regions of galaxies due to tidal debris from the accretion and disruption of 
an expected population of ~ 100 — 200 luminous satellites. Gauthier et al. (2006) and 
Kazantzidis et al. (in prep) further show how the accretion of such a population will heat 
and restructure the thin disk, producing additional "low latitude" debris. I summarize 
here our current understanding of the substructure seen around M31: 

• Deep HST/ACS CMDs reaching well below the horizontal branch have been obtained 
for 8 regions in the far outskirts of M31, including seven regions of visible substructure. 
Many of the CMDs exhibit different morphologies, suggesting variations in the mean age 
and metallicity of the constituent stars (Ferguson et al. 2005). In all cases, the metallicity 
inferred from the RGB colour is significantly higher ([Fe/H]> —0.7) than that of typical 
present-day low mass Local Group dwarf satellites ([Fe/H]< —1.5), implying that such 
objects are unlikely to be the progenitors of the observed substructure. 

• All substructure fields lying near the major axis contain stars of age < 1.5 Gyr 
(Faria et al 2007; Richardson et al., in preparation); this includes the NE clump field, 
which lies at a projected radius of 40 kpc. Additionally, the major axis substructure fields 
are dominated by a strong rotational signature, similar to the HI disk, with a modest 
velocity dispersion (Ibata et al. 2005). Taken together, these results suggest that the 
"low-latitude" structure in M31 results primarily from perturbed outer thin disk. 

• The combination of line-of-sight distances and radial velocities for stars at various 
locations along the giant stream constrains the progenitor orbit (e.g. Ibata et al. 2004; 
Fardal et al. 2006). Currently- favoured orbits do not easily connect the more luminous 
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inner satellites (e.g. M32, NGC 205) to the stream however this finding leaves some 
remarkable coincidences (e.g. the projected alignment on the sky, similar metallicities) 
yet unexplained. The giant stream is linked to another overdensity, the diffuse feature 
lying north-east of M31's centre, on the basis of nearly identical CMD morphologies and 
RGB luminosity functions; orbit calculations suggest this connection is indeed likely. 

3. Diffuse Stellar Halos 

Tidal debris from accretion events will disperse over time and streams from ancient 
events will have long merged to produce a diffuse stellar halo (e.g. Abadi et al. 2005; 
Bullock & Johnston et al 2005). If hierarchical growth is the main mode of mass assembly 
then such stellar halos should be a generic feature of galaxies. Unfortunately, there are 
very few observational constraints on the nature and ubiquity of stellar halos around 
galaxies (see however Zibetti et al. 2004). Even our understanding of the M31 and M33 
stellar halos was extremely poor until recently. I summarize some key results here: 

3.1. MSI 

Using data from the INT/WFC survey, Irwin ct al. (2005) combined diffuse light surface 
photometry with resolved star counts to probe the minor axis profile of M31 to a radius 
of ~ 55kpc (effective fiy ~ 32 mag/D"). The profile shows an unexpected flattening (rel- 
ative to the inner R 1 / 4 decline) at a radius of ~ 30kpc, beyond which it can be described 
by shallow power-law (index « —2.3), possibly extending out to 150 kpc (Kalirai et al. 
2006). This compares favourably with the Milky Way halo, which exhibits a power-law 
index of —3.1 in volume density (Vivas et al. 2006). The discovery of a power-law com- 
ponent which dominates the light at large radii in M31 has profound implications for the 
interpretation of all prior studies of the M31 "stellar halo" . Since these studies generally 
targetted regions lying within 30 kpc along the minor axis, they most likely probed the 
extended disk/bulge region of the galaxy and not the true stellar halo. 

Kcck/DEIMOS spectroscopy has been used to study the kinematics and metallicities of 
stars in the far outer regions of M31. By windowing out the stars which corotate with the 
HI disk, Chapman et al. (2006) have detected an underlying metal-poor ([Fe/H]~ —1.4), 
pressure-supported (er ~ 100 km/s) component (see also Kalirai et al. 2006). Although 
it has yet to be proven that this component is the same one which dominates the power- 
law profile at very large radius, the evidence is highly suggestive. Despite many years 
of thinking otherwise, it thus appears that M31 does indeed have a stellar halo which 
resembles that of the Milky Way in terms of structure, metallicity and kinematics. 

3.2. M33 

Studies of the M33 stellar halo have also had a somewhat checkered history. Mould & 
Kristian (1986) measured [M/H]^ —2.2 in a field located at 7 kpc along the minor axis 
of M33 and, for the better part of two decades, this metallicity was generally assumed 
to reflect that of the M33 halo. Recent work has found a significantly higher metallicity 
for stars in this same field and, at the same time, suggested that the field is actually 
dominated by the outer disk and not the stellar halo (e.g. Ticdc ct al. 2004, Ferguson et al. 
2007 in prep). The detection of a power-law structural component in the outskirts of M33 
has so far proved elusive, although the RGB clearly becomes narrower and more metal- 
poor in these parts (Ferguson et al. 2007 in prep). A recent Keck/DEIMOS spectroscopic 
study has targetted two fields in the outskirts of the galaxy, located at ~ 9 kpc along the 
major axis (McConnachie et al. 2006). Although the dominant kinematic component in 
these fields exhibits strong rotation, there is tentative evidence for an additional low-level 
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metal-poor ([Fe/H]~ —1.5) component centered at the systemic velocity. While further 
work is clearly needed to elucidate the nature of this structure, the estimated velocity 
dispersion (cr ~ 50 km/s) supports an association with M33's true stellar halo. 

4. Future Outlook 

Quantitative study of the faint outskirts of galaxies provides important insight into the 
galaxy assembly process. Much recent work has focused on our nearest large neighbours, 
M31 and M33. While both appear to show evidence for a metal-poor, pressure-supported 
stellar halo (similar to that of the Milky Way), only M31 shows evidence for recent accre- 
tion. In order to put these results in context, it is necessary to establish the properties of 
a larger sample of galaxies, spanning a range in both host galaxy luminosity and Hubble 
type. With 8-m class telescopes, it is possible to map the RGB populations in the out- 
skirts of galaxies to distances of < 5 Mpc; with HST, this work can feasibly be extended 
to > 10 Mpc. Spectroscopic characterization of resolved RGB populations beyond the 
Local Group is far harder and must await the arrival of 30-m class telescopes. 
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